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We address local quantum estimation of bilinear Hamiltonians probed by Gaussian states. We evaluate the rel- 
evant quantum Fisher information (QFI) and derive the ultimate bound on precision. Upon maximizing the QFI 
we found that single- and two-mode squeezed vacuum represent an optimal and universal class of probe states, 
achieving the so-called Heisenberg limit to precision in terms of the overall energy of the probe. We explicitly 
obtain the optimal observable based on the symmetric logarithmic derivative and also found that homodyne 
detection assisted by Bayesian analysis may achieve estimation of squeezing with near-optimal sensitivity in 
any working regime. Besides, by comparison of our results with those coming from global optimization of the 
measurement we found that Gaussian states are effective resources, which allow to achieve the ultimate bound 
on precision imposed by quantum mechanics using measurement schemes feasible with current technology. 

PACS numbers: 03.65.Ta, 42.50.Dv 
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I. INTRODUCTION 

In this paper we address quantum estimation of unitary op- 
erations for continuous variable systems. In particular we an- 
alyze the estimation of the interaction parameter 6 for uni- 
taries of the form Ug = exp{— i9G} where G is a linear or 
bilinear bosonic Hamiltonian of the form G = a^b + aw, 
G = at&t + ab, or G = a t2 + a 2 , [a, a*] = 1 and [b, &+] = 1 
being mode operators. We are interested in evaluating the ul- 
timate bound on precision (sensitivity), i.e the smallest value 
of the parameter that can be discriminated, and to determine 
the optimal measurement achieving those bounds. 

As a matter of fact, linear and bilinear interactions for 
bosonic systems are a key ingredient for the development of 
continuous variable quantum information processing [2, IH 
4[. They are usually realized by means of parametric pro- 
cesses, as single- and two-mode squeezing, or by linear op- 
tical elements such as phase-shifting and two-mode mixing. 
The precise characterization of linear optical gates is also of 
interest in interferometry S 0. §], absorption measurement 
Jit] and characterization of detectors ifloll . 

In general, interaction parameters cannot be directly ac- 
cessed experimentally, and the estimation process consists in 
probing the interaction by a known quantum signal qq, which 
is measured after the interaction (see Fig. [TJ. The relevant 
constraint in the optimization of those schemes concerns the 
total energy of the probe, which should be kept as low as 
possible to avoid any possible modification or degradation of 
the gate itself. Overall, the problem we are facing is that of 
devising the optimal measurement, i.e. a positive operator- 
valued measure (POVM) {E x } x& x, to be performed on the 
probe Qg — UgQoUl after the interaction, at fixed energy 
N = Tr[f?o J2j n j] °f the incoming signal, J^j n j being the 
total number operator of the involved modes. 

The above problem may be properly addressed in the 
framework of quantum estimation theory (QET) lTTTl[l2i[T3ll . 
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FIG. 1 : General scheme for the indirect estimation of the the unitary 
Ue probed by the signal go. 



which provides analytical tools to find the optimal measure- 
ment according to some given criterion. In turn, there are two 
main paradigms in QET: Global QET looks for the POVM 
minimizing a suitable cost functional, averaged over all pos- 
sible values of the parameter to be estimated. The result of 
a global optimization is thus a single POVM, independent on 
the value of the parameter. On the other hand, local QET 
looks for the POVM maximizing the Fisher information, thus 
minimizing the variance of the estimator, at a fixed value of 
the parameter lfl4lll5ll . Roughly speaking, one may expect lo- 
cal QET to provide better performances since the optimization 
concerns a specific value of the parameter, with some adap- 
tive or feedback mechanism assuring the achievability of the 
ultimate bound [16]. Global QET has been mostly applied to 
find optimal measurements and to evaluate lower bounds on 
precision for the estimation of parameters imposed by unitary 
transformations. For bosonic systems these include single- 
mode phase ifnl [l8ll . displacement [19], squeezing I20I l2ltl 
as well as two-mode transformations, e.g. bilinear coupling 
[9J]. Local QET has been applied to the estimation of quan- 
tum phase [24] and to estimation problems with open quan- 
tum systems and non unitary processes i25ll : to finite dimen- 
sional systems 12611 . to optimally estimate the noise parameter 
of depolarizing [27] or amplitude-damping lf28ll . and for con- 
tinuous variable systems to estimate the loss parameter of a 
quantum channel 112911 . 

In this paper we consider the estimation the interaction pa- 
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rameters of bilinear bosonic Hamiltonians from the perspec- 
tive of local QET. In particular, we focus our attention to mea- 
surement schemes as in Fig. Q]with the probe state chosen 
within the set of Gaussian states 01 0, HHH H3l, which rep- 
resents a class of signals achievable with current technology. 
We evaluate the relevant quantum Fisher information (QFI) 
and derive the ultimate bound on precision. Upon maximiz- 
ing the QFI we found that single- and two-mode squeezed 
vacuum represents an optimal and universal class of probe 
states, achieving the so-called Heisenberg limit to precision 
in terms of the overall energy of the probe. Remarkably, by 
comparison with results coming from global optimization of 
the measurement |0, |2(| [2j]] we found that Gaussian states 
are effective resources, which allow to achieve the ultimate 
bound on precision. Besides, we found that homodyne detec- 
tion assisted by Bayesian analysis may achieve near-optimal 
sensitivity in any working regime. 

The paper is structured as follows: in the next Section we 
briefly review local quantum estimation theory with some re- 
marks on the implementation of the optimal measurements. In 
Section|III]we evaluate the optimal measurements and the cor- 
responding bounds on precision for the local estimation of bi- 
linear couplings using Gaussian probes. In Section[IV]we ad- 
dress estimation of squeezing using homodyne detection and 
Bayesian analysis and show that near-optimal precision may 
be achieved in any working regime. In Section [V] we com- 
pare our results with those coming from global estimation and 
close the paper with some concluding remarks. 



II. LOCAL QUANTUM ESTIMATION THEORY 

In this section we review some concepts of local quantum 
estimation theory Iq3l l23tl which will be used in the rest of 
the paper. As a matter of fact, many quantities of interest in 
different branches of physics cannot be directly accessed ex- 
perimentally, either in principle, as in the case of field mea- 
surement 13 in . or due to experimental impediments. In these 
cases, one has to indirectly estimate the value of those phys- 
ical parameters by measuring a different observable, some- 
how related to the quantity of interest. This indirect procedure 
of parameter estimation implies an additional uncertainty for 
the measured value, that cannot be avoided even in optimal 
conditions. The aim of quantum estimation theory is to op- 
timize the inference procedure by minimizing this additional 
uncertainty. In the classical theory of parameter estimation 
the Cramer-Rao Bound i32ll establishes a lower bound for the 
variance of any unbiased estimator 9 of the parameter 9. This 
lower bound is given by the inverse of the so-called Fisher 
Information (FI): 



A9 2 > 



1 



F{9) 



where the Fisher Information is defined as 

' d\np(x\> 



F{9) =Y,P( X \°) 
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(1) 

(2) 
(3) 



Here 9 is the parameter to be estimated, and x denotes the 
outcome of the measurement of the quantity X related to 9. 
The notation p(x\9) indicates the conditional probability of 
obtaining the value x when the parameter has the value 9. 

A quantum analogue to Eq. (01 may be found starting from 
the Born rule 



p(x\9) = tr[E xPe ] 



(4) 



where E x are the elements of a positive operator-valued mea- 
sure (POVM) and pg is the density operator, parametrized by 
the quantity of interest, describing the quantum state of the 
measured system. The Fisher Information is then rewritten as 



m = E 



ReTi[p e E x K 6 



(5) 



where we introduced the Symmetric Logarithmic Derivative 
(SLD) Kg, which is the self-adjoint operator defined as 



dpe 
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(6) 



It can then be shown flilfTsll that the Fisher Information (0 is 
upper bounded by the so-called Quantum Fisher Information 
(QFI): 



F < H = Ti[peA 2 e } 



(7) 



In turn, the quantity 1 /H represents an ultimate lower bound 
on precision for any quantum measurement (followed by any 
classical data processing) aimed to estimate the parameter 9. 
The SLD is itself an optimal measurement, that is, using the 
POVM E x obtained from the projectors over the eigenbasis of 
Kg we saturate the inequality (|7]). 

In this work we will focus on systems where the depen- 
dence of pg from the parameter 9 is generated by a fam- 
ily of unitary transformations: pg = UgpoU\ where Ug — 
exp(— i#G), G is the Hamiltonian that generates the trans- 
formation and go is a given quantum state used to probe the 
Hamiltonian process. In this case it is possible to obtain an 
explicit formula for the SLD operator and the QFI. At first we 
take the eigenbasis of p Q : p = I^.PA#fc)(^fe| From © we 
can rewrite Kg in this basis as follows 



Kg = 2i E G 



jk 



jk 



Pj ~ Pk 
Pj, + Pk 



Ug\^ 3 )(MUl 



(8) 



where Gjk = (ipj \G\ipk) are the matrix elements of the gen- 
erator G. Eq.® shows that Kg depends on 9 only through 
the unitary transformation Ug. As a consequence it is possi- 
ble to define the operator Ao, independent from 9, such that 
= UgK Q Ug. It also follows that the quantum Fisher in- 
formation is independent from 9. In fact, H = Tr[peAg] = 

Tv[UepoUlUgKlul] = Tr[p A§]. Explicit formulas to calcu- 
late H may be given in the eigenbasis of po 



i ik 



Pn - Pk 
l P n +Pk 



C 2 
^ n k 



4(G 2 



n k 



PkPn 
Pn + Pk 



G n kGk 



(9) 
(10) 
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As we will see in the following, situations with a probe de- 
scribed by a pure state p Q = \ipo)('4>o\ are of particular in- 
terest. In those cases the QFI reduces to the variance of the 
generating Hamiltonian G, i.e. H = 4 AG 2 . In addition, for a 
pure state we have p 2 e — pg and thus Aq = 2i[po, G] i.e. 



An = 2 



J2( G ^o)(ipk\-G ko \^ k ){il>o\ 



(11) 



III. ESTIMATION OF BILINEAR COUPLINGS 

In this Section we address the case of local estimation of 
various bilinear couplings (single- and two-mode squeezing, 
two-mode mixing) using Gaussian probes at fixed energy. 



A. Single-mode squeezing 

Here we consider the estimation of the parameter 9 im- 
posed by the unitary transformation exp (— iOG), where G is 
the generating Hamiltonian 



G 



1 



(a t 2 + a 2 ) 



(12) 



We analyze the precision achievable in the estimation of 9 by 
using different classes of (Gaussian) probe states. The mea- 
surement aimed to estimate 9 is made on the transformed state 



pe = exp (—i9G) po exp (—i9G) 



(13) 



At first we analyze the case of a Gaussian pure probe i.e 
a squeezed coherent state of the form p = IV'oXV'ol with 
\tpo) > = S(r)D(a)\0), where 



D(a) 
S(r) 



exp \aa 



exp 



t 



a a 



(14) 
(15) 



and where, without loss of generality, we have chosen a real 
squeezing parameter r and a complex displacement a = xe 1 ^ . 
Since po is a pure state, the QFI will be given by 



H = 4AG 



((a 



t2 



,t2 



2\2 



Upon evaluating all the expectation values we obtain: 



AG 



-x 2 cos 2cj> sinh 2r + (2N + 1) sinb 2 r + N- 



(16) 



2 

(17) 



where N = (aJa) = x 2 + sinh 2 r denotes the overall energy 
of the probe signal. The signal optimization corresponds to 
the maximization of H over the state parameter with the con- 
straint of fixed N . The phase is a free parameter since it 
does not influence the total energy. The choice cos 2<f> = — 1 
maximizes H leading to 

H = 4(A - sinh 2 r) sinh 2r + 4(2A + 1) sinh 2 r + 4A + 2 

(18) 



which grows monotonically with sinh r and achieve its max- 
imum 



8N 2 



8N - 



(19) 



for sinh r — N and a — 0, corresponding to a squeezed 
vacuum probe. Thus, to obtain the maximum accuracy in the 
estimation of 9 it is more efficient to use all the energy in 
squeezing rather than field amplitude. 

In order to see the effects of mixing we have also considered 
a class of probes made by squeezed thermal states 



1 



Po 



S(z)\k)(k\SHz) 



(20) 



where the squeezing z — re 1 ^ is a complex number. We are 
now dealing with a mixed state; the corresponding QFI is thus 
given by The state vectors of the diagonal basis of po and 
their associated probabilities are 



S(i9)\k) 
1 



Pk 



n + 1 \n + 1, 
The matrix elements of the generator G are 



(21) 
(22) 



G jk =(k\SHz) 

_ 1 r 
~ 2 



,.t2 



S(z)\k) 



V(j + l)(j + 2)( Ai 2 + v* 2 ) 5, 



j+2,k 



+ ^/(k + l)(k + 2)(p 2 + v 2 ) Sj, k+2 
+ {2k + l)p{v + v*)8 :jM 



(23) 



where p — cosh r, v — e 1 ^ sinh r. From this and (0 we get 



H =2 (cosh 4 



sinh 4 r + 2 cos 2cf> sinh 2 r cosh 2 r) 



4n 2 + 4n + 1 
2n 2 + 2n + 1 
The energy constraint is now given by 



N 



+ (2n + 1) sinh 2 r 



Maximization over the free parameter <j) leads to <p - 
turn to 

(4n 2 + 4n+ 1)(4A 2 + 4A+ 1) 



H 



(2n 2 + 2n+l)(2n+l) 2 



(24) 

(25) 
and in 

(26) 



The maximum of this function is found when n = 0: again 
we are led to squeezed vacuum. 

As we have already discussed, the optimal measurement, 
i.e. when the Fisher Information is equal to the QFI, is realized 
by the SLD A. For squeezed vacuum probes we may use Eq. 
(fTTT l and obtain 



A = i\/2(2N + 1 



)ff(r){|0)<2|-|2)(0|}st(r) 



(27) 



Summarizing, the most convenient way of estimating a 
squeezing parameter is to probe the transformation by a 
squeezed vacuum probe. The corresponding QFI scales as 
H ~ 8A 2 in terms of the overall energy of the probe. 
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B. Two-mode mixing 

Here we consider the case where the generator G is the two- 
mode mixing Hamiltonian: 



G = a?b + atf 



(28) 



Let us first consider a probe state made by factorized squeezed 
thermal states: 



Pa = [S a {r) ® S b {s)]v a ® i*[Sj(r) ® s£(s)] 
where v ajb are the density matrices of thermal states: 



1 



Vk = 



(n fe + 1) 



n fe + l 



|n)(n| 



(29) 



(30) 



For a two-mode system the formula (O for the Quantum 
Fisher Information becomes 



H = 4 Pi 



jkmn 



Pjk Pmr< 
Pjk H~ Pmr, 



G 



jkmnGm nj k 



(31) 



where pkn — PkPn, the thermal coefficients d22l . The Heisen- 
berg evolution of the mode operators 



Sl(r)Sl(s) (a^b + atf) S a (r)S b {s) = cosh(r + s)(a t 6 + a6 f ) + mnh(r + s)(ab + aW) 
allows to calculate the matrix elements of G 

G jta =0-,fc|S+(r)5 6 t (s) (a^b + atf) S a (r)S b (s)\ m, n) 

= cosh(r + s) (j (m +l)(k + l)5j= m +i^n=fc+i + + l)(n + l)5 m =i+i^fc=n+i) 
+ sinh(r + s) (V (i + + l)(5 m =j+i^=fc+i + \/(m + l)(n + l)^ =m+ i(5fc=n+i) 
The resulting QFI reads as follows 

(n a -n b ) 2 (n a +n b + l) 2 \ (n a - n b f 



H = 4 



sinh (r + s) 



2n a n b + n a + n b 2n a n b + n a + n b + 1 / 2n a n b + n a + n b 
I 



(32) 



(33) 



(34) 



The total photon number of the system is given by the sum 

N = n a + n b + (2n a + 1) sinh 2 r + (2n b + 1) sinh 2 s 

(35) 

The QFI d34i > has no point of gradient zero that is compatible 
with the energy bound d35t . Since it is a continuous function, 
to find its maximum we need to investigate its value at the bor- 
ders of its domain. Let us first consider the case n a — n b = 0, 
i.e. a probe made by two disentangled squeezed vacuums. The 
energy and the QFI become respectively 

(36) 
(37) 
s, which 

(38) 



whose maximum is 



N 



sinh 2 r + sinh 2 s 



H = 4 sinh (r + s) 

The maximum of this function is reached when r 

gives 

Hi = 47V 2 + 87V 



The second possible case is given by two thermal states, when 
r = s = 0. The QFI becomes 

4(JV- 2n b f 



H = 



N + 2(N - n b )n b 



(39) 



H2 = 47V when n a = or n b = 



(40) 



i.e. when one of the states is at zero temperature. The last pos- 
sible combination is given by a thermal state and a squeezed 
vacuum, for r = 0, n b — 0. Energy and QFI reduce to 

7V = n a + sinh 2 s (41) 

H = 4 [(2n„ + 1) sinh 2 s + n a ] = 4[7V + 2n a (N - n a )\ 

(42) 

The optimal QFI is obtained when the energy is equally dis- 
tributed between the thermal state and the squeezed state, 

• i 2 N 

n a = sinn s = -5-: 



H* = 27V 2 



■47V 



(43) 



Thus we see that the maximum Fisher information is obtained 
using two equally squeezed vacuums. Since this is the combi- 
nation of two pure states, we can use ([TO to obtain the SLD 
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that realizes the optimal measurement: 

A = 2i^N(N + 2)5 a 5 6 (|0,0)(l,l| - |1, 1) (0, O^S* 

(44) 

In order to investigate the role of entanglement in the estima- 
tion procedure we consider the probe prepared the state 

Po = l^ooXlM (45) 

l^-fc) = A)» = U(cj>)T(X)\j,k) (46) 

where 

U{<t>) = exp[-^(a6 t + a f 6)] (47) 

T(A) = exp[-iA(a6 + aW)] (48) 

The probe is transformed into pg = e~ l9G p^e l9G where again 
we are using the generator (l28l l. Since we are dealing with a 
pure state, the QFI is 

H = 4AG* 2 = 16 cosh 2 |A| sinh 2 |A| (l - 4 cos 2 sin 2 <f>) 

(49) 



where we have used the Heisenberg evolution of the mode 
operators. The energy constraint is given by 



N = (a f a) + (6+6) = 2 sinh 2 |A| 



thus the QFI can be rewritten as 



H = (4iV 2 + 8N) (1 - 4 cos 2 sin 2 



(50) 



(51) 



Since d^N = 0, we can freely choose a value for cf>, in order 
to maximize H . The maximum Fisher information is obtained 
for cos 40 = 1 and corresponds to H = 4N 2 + 8N, i.e. no 
improvement is obtained using an entangled probe. The SLD 
operator that realizes the optimal measurement is found using 
CD: 



A = 2i^2iV(iV+l){|V>oo>(V'2o| + iV-ooXlM - |Vao)<Vto| - I V^02> (V^oo I } 



(52) 



r 



C. Two-mode squeezing 

The procedure used for the case of two-mode mixing may 
be analogously applied when the generator G is given by the 
two-mode squeezing Hamiltonian: 



G = ab + aW 



(53) 



First we analyze the case of an initial density matrix, see I 
that describes two disentangled squeezed thermal states. The 
same steps done to obtain d34] l can be repeated, using the 
Hamiltonian d53l l instead of ( 1281 . The QFI for this particu- 
lar case is thus given by 



H = 4 



sinh 2 (r + s) 



(ni - n 2 ) 2 



(ni + n 2 + 1) 



(ni +n 2 + l) 2 



1n\n 2 + n\ + ri2 1ri\n 2 + n\ + n% + 1 / 2n\n 2 + n\ + n 2 + 1 
I 



(54) 



The maximum of this function is once again obtained when The same can be done for the case of a probe such as d45l ). 
fix = n 2 = and r = s, i.e. when the probe is made by two The corresponding QFI is given by 



equally squeezed vacuum states. This max is 

# max = 4(2iV+l) 2 



(55) 



H =8 cosh 2 | A| [(cos 2 - sin 2 (ff cos(2 arg A) sinh 2 | A| 
+2 sinh 2 |A| + 1] 



(57) 

= 4N 2 + %N is 

where N = 2 sinh r. The corresponding SLD reads as follows achieved when cos(arg2A) = 1 and cos 20 = 1 and using 

the SLD 



The maximum Fisher Information H„ 



A = 2i(N+l)S a S b (\0,Q)(l,l\ - |1, 1X0,01)^ (56) A = 2i(2JV+l)(|V>oo>fcM - \^n)(^oo\] 



(58) 
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IV. ESTIMATION OF SQUEEZING BY HOMODYNE 
DETECTION 

In Section |III] we have shown that squeezed vacuum is the 
optimal reference Gaussian state to estimate the parameter 
of a squeezing transformation. However, the optimal mea- 
surement maximizing the QFI, that is the SLD, is not realiz- 
able with current technology. It is thus of interest to inves- 
tigate whether a feasible measure may be used to effectively 
probe the perturbed squeezed vacuum. We focus to the case of 
single-mode squeezing estimation; an analogue analysis may 
be performed for two-mode operations. Our approach is to 
exploit homodyne detection to measure field-quadrature: 



1 



[ae +a'e J 



(59) 



and inferring the squeezing parameter through the results ob- 
tained with multiple homodyne measurements. The homo- 
dyne probability p(x\9) is given by 



p(x\6) = Tr[p e 11,(0)] 



(60) 



n x = \x)eg(x\ being the spectral measure of the quadra- 
ture t59) . The resulting distribution for a squeezed vacuum 
to which an unknown squeezing has been applied, is a zero 
mean (Tr [pg x a ] — 0) Gaussian distribution 



p(x\9) 



ex p(-^2! > (61) 



with variance (see the Appendix for details on the derivation) 



Y? g = cos{2u)^N{N + 1) 

N + [cosh(20) + sin(2a) sinh(20)] (62) 

The reason to choose homodyne detection is that the classical 
Fisher information (f3]) of the homodyne distribution p a (x\9) 
may be optimized over a in order to achieve the same scal- 
ing as the QFI versus the energy of the probe. Being a\ = 

arg max Q F a (0) we have 



cos ai 



F ai (e) N ^8Ni 



y/NW+T) 



(l + 2iV)cosh0-sinh0 



(63) 
(64) 



This means that homodyne detection with optimized phase a 
is a good candidate to achieve ultimate bounds to precision, 
as far as it saturates the classical Cramer-Rao bound. Indeed, 
Von Mises-Bernstein-Laplace theorem ensures that Bayesian 
a posteriori distribution p(9\{x}m), representing the proba- 
bility of the squeezing to be given the homodyne sample 
{x}m, converges asymptotically to a Gaussian distribution, 
centered in the true value with variance saturating the Cramer- 
Rao bound. In other words, Bayesian estimators are asymptot- 
ically unbiased and efficient. In the following, we thus discuss 



in some details estimation of squeezing by homodyne detec- 
tion and Bayesian analysis. We consider a large number M 
of homodyne measurements on repeated preparations of the 
same system. Since the measurements are independent, the 
a posteriori distribution is proportional to the product of the 
single data distribution 



M 



p{9\{x} M )^\{p(9\x k ) = Yl 



k=l 



k=l 



P(x k \9)p(9) 
p{xk) 



(65) 



where we repeatedly used the Bayes Theorem. p(9) is the a 
priori distribution of the parameter, p(x) the overall probabil- 
ity of the outcome x, while p(x\9) is the probability to obtain 
the outcome x when the squeezing parameter is 0. The proba- 
bility p(9\{x}j\f) has to be normalized, Eq.d65b thus rewrites 
as 

p(e\{x} M ) = ^p(erH^f^ (66) 

A £A p{xk) 

where A is the normalization constant given by 



A = 



-f-OC 



W TT P( x k\Q) 



p(o) m n 



k=l 



p(x k ) 



(67) 



We assume to have no a priori information on the squeez- 
ing i.e. we take p(9) as a uniform function. Notice also 
that the product of the distributions p(x k ) does not depend on 
and it cancels out due to normalization. Finally, since we 
wish to perform a large number M ^> 1 of measurements, 
the product in d66*l l will contain many repeated elements: each 
outcome x is obtained a number of times proportional to its 
probability p(x\9*), being 0* the true (and unknown) value of 
the squeezing parameter. We can then re-order the product so 
that its index now runs through all possible values of x: 



\Mp{x\8") 



p(9\{x} M )^jY[p(x\9Y 

X 

= -^exp |m J p(x\9*)lnp{x\9)dx^ (68) 

where we have taken a limit to the continuum for the variable 
x. The integral in (|68T > can be solved leading to 



p(x\9*) hip(x\9)dx 



§ + ln(27r£2) 



(69) 



where we have introduced the short notation = S^, . Over- 
all, we obtain 



P{0\{x}u) = j 



-M/2 



(70) 



where we have redefined A so to include all terms independent 
from 0. The mean of the a posteriori distributionp(0|{a;} m) 
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is our estimator and the variance A9 2 the corresponding con- 
fidence interval 



f M6p(0\{x} M ) 

' — oo 
p+ca 

A9 2 = ' JDIa m2 

J — OO 



dd(9-e) 2 P (0\{x} M ) 



(71) 



(72) 



An optimal value for the homodyne phase a is obtained upon 
minimizing the variance of the a posteriori distribution. Be- 
sides the value a% reported above we found that optimal scal- 
ing (oc A/ -1 A~ 2 ) of the variance may be achieved also for 
the phase value 



OL2 



-sign( 



scch(26>*)sinh 2 9* 



which, remarkably, is independent on the probe energy A (in- 
deed, we have a% = a-i + 0(l/N)). 
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FIG. 2: (Color online) Left: Rescaled a posteriori distribution qu (9) 
for M = 5, N = 40 (black), M = 10, N = 20 (blue), M = 
20, N = 10 (red), M = 40, N = 5 (green). Right: LogPlot of the 
rescaled a posteriori distribution, for the same values of the parame- 
ters. 

In Fig. |2] we report the rescaled distribution qui (9) = 
p(9\{x}m)/(MN) for different values of the probe en- 
ergy and the number of measurements, we also report 
p(9\{x}m)I (MN) in a logarithmic scale to enlighten the dif- 
ferences in the distribution tails. As it is apparent from the 
plots the relevant parameter is the energy of the probe. For 
highly excited probes, i.e for JV > 1, we expand I] 2 , as 



A 



1 
2 

cos(2o! 
8A 



[cos(2a) 



O 



1 

A 2 



cosh(26») + sin(2a) sinh(20)] 
(73) 



and neglect all orders scaling as A 2 or higher. Upon choos- 
ing the homodyne phase &2 we have 



a _ scch(2g*) 
8A 

E 2 , ~ sech(20*) 



(2A+l)sinh 2 (<9-6>*) + 



8A 



(74) 
(75) 



Upon substituting ( f74b and (Tt3T > into ( l70b we see explicitly that 
p(9\{x}m) — p{9— 9*\{x}m) and that the estimator is indeed 
unbiased, i.e 9 = 9*. We also found that the variance is inde- 
pendent from the true value of the squeezing 9*: Numerical 

computation shows that the variance A9 2 scales as ~ - 

F 4Af A^ 2 

for large N, that is, apart from a factor two, the same scal- 
ing of the inverse of the QFI ( fl"9l l. Notice that the optimal 
phase «2> depends on 9* , which is the unknown parameter 
that we are trying to estimate. This is consistent with the lo- 
cal nature of the estimator procedure. From a practical point 
of view this means that some kind of feedback mechanism or 
adaptive technique should be employed to adjust the phase of 
the homodyne detector 11161 13311 . We conclude that homodyne 
detection with Bayesian analysis is a robust and accurate es- 
timation technique for the squeezing parameter. Remarkably, 
this scheme may be implemented with current technology. 



V. CONCLUSIONS 

In this paper we have addressed local quantum estimation 
of bilinear Hamiltonians probed by Gaussian states. We eval- 
uated the relevant quantum Fisher information (QFI) thus ob- 
taining the ultimate bound on precision. Upon maximizing 
the QFI we found that single- and two-mode squeezed vac- 
uum represent an optimal and universal class of probe states, 
achieving the so-called Heisenberg limit to precision in terms 
of the overall energy of the probe. For two-mode operations 
no improvement may be obtained using entangled probes. 

It is worth noting that the Heisenberg scaling A9 ~ A -1 
in terms of the overall energy of the probe may be achieved 
also using global quantum estimation techniques (see e.g. 
for the case of two-mode mixing). In that case, however, op- 
timization of the probe have been performed over the whole 
set of quantum states, not focusing on Gaussian states. In 
turn, this means that Gaussian states are effective resources, 
which allow to achieve the ultimate bound on precision im- 
posed by quantum mechanics using measurement schemes 
feasible with current technology. This has been confirmed by 
a Bayesian analysis applied to the estimation of squeezing by 
homodyne detection, which achieves near-optimal sensitivity 
in any working regime, i.e for any (true) value of the squeez- 
ing parameter. For the estimation of squeezing, Heisenberg 
scaling for Gaussian probes has been also found exploiting 
global strategies 12111 . In that case, however, though the mea- 
surement does not depend on the value of the parameter, there 
is a strong dependence on the probe states. We have also ex- 
plicitly obtained the optimal observables based on the sym- 
metric logarithmic derivative, which however do not corre- 
spond, in general, to a feasible detection scheme. 

We conclude that Gaussian states and Gaussian measure- 
ments assisted by Bayesian analysis represent robust and ac- 
curate resources for the estimation of unitary operations of 
interest in continuous variable quantum information. 



J 



VI. APPENDIX 

Here we show how Eq.(l62b is obtained. We start from the identity 

S' (z)aS(z) — (ia + va' 

where /i = cosh | z | and v — e % arg z sinh | z | . In turn this leads to 

{r)S^ (i9)aS(i9)S(r) = (acoshr + sinhr) cosh (9 + i{a i coshr + asinhr) sinh 9 

= a (cosh r cosh 9 + i sinh r sinh 9) + (sinh r cosh 9 + i cosh r sinh 9) 

and then 

+ /i.e.} 
-/i.e.} 



&(r)&(iff)x a S(i8)S(r) = i{ 

-5M< 



( ch r ch 6* — i sh r sh 0) + a( sh r ch 6* — z ch r sh ( 



' ch r ch 9 — i sh r sh ( 



(shrch6* + zchrsh.0)e 



(76) 



(77) 



(78) 



When the square of this operator is averaged in the vacuum (0|...|0), only one of the four terms a 2 , a^a, aa^ and aJ 2 survives, 
namely (0|<z<r |0) = 1. The equation then simplifies to 

S 2 , =ije 2M (chrch6l-ishrsh60(shrch(9- ichrsh9) + e" 2M (shrch6» + ichrsh6>)(chrch6» + ishrsh6») 
+ ( ch r ch 9 — i sh r sh 9) ( ch r ch 9 + i sh r sh 9) + ( sh r ch 6> — i ch r sh 9) ( sh r ch 9 + i ch r sh 9) j 
= i| sinh(2r) cos(2a) + cosh(2r)[cosh(2<9) + sin(2a) sinh(2^)]| 

= cos(2a) T/N(NTTj + (n + ^ [cosh(26») + sin(2a) sinh(20)] (79) 



where we used N = sinh r. 
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